A long gamma-ray burst from a merger of compact objects 
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Gamma-ray bursts (GRBs) are flashes of high-energy radiation arising from energetic 
cosmic explosions. Bursts of long (>2 s) duration are produced by the core-collapse of massive 
stars!, those of short (<2 s) duration by the merger of two neutron stars? (NSs). A third class 
of events with hybrid high-energy properties was identified*, but never conclusively linked 
to a stellar progenitor. The lack of bright supernovae rules out typical core-collapse 
explosions*>°, but their distance scales prevent sensitive searches for direct signatures of a 
progenitor system. Only tentative evidence for a kilonova has been presented’. Here we 
report observations of the exceptionally bright GRB211211A that classify it as a hybrid event 
and constrain its distance scale to only 346 Mpc. Our measurements indicate that its lower- 
energy (from ultraviolet to near-infrared) counterpart is powered by a luminous (~10*” erg 
s') kilonova formed in the ejecta of a compact binary merger. Our work demonstrates that 
compact binary mergers can power high-energy transients of over a minute long duration, 
and that these represent a new, albeit rare, electromagnetic counterpart of gravitational 


wave radiation. 


On 11 December 2021 at 13:59:09 Universal Time (UT; hereafter To), NASA’s Neil Gehrels Swift 
observatory (hereafter Swift) discovered GRB211211A® as an extremely bright burst with a 
duration of over 50 s (Extended Data Figure 1). The burst was independently observed by the 
Fermi, INTEGRAL, and CALET satellites. Its optical, ultraviolet (UV) and X-ray counterparts were 
localized within minutes, close to a nearby galaxy, SDSS J140910.47+275320.8 (G1 in Figure 1), 
at a distance of 346 Mpc (Methods). Spectroscopic observations of the optical counterpart showed 
a featureless continuum? and did not allow for a direct measurement of the GRB distance scale. 


However, when combined with the detection of a bright UV counterpart and the lack of any 


underlying host galaxy down to F/60W> 27.6 AB mag, these observations point to a low-redshift 
origin for GRB211211A (z<1.5 at the 99.9% confidence level; Methods). The most likely birthsite 
is in the outskirts of the nearby galaxy G1, at a projected physical offset of 8.00 +/- 0.04 kpc from 
the galaxy’s nucleus. This association is also supported a) by probabilistic arguments, being the 
chance alignment between the GRB and the bright G1 galaxy only 1.4%, b) by the uncommon 


brightness of the prompt gamma-ray emission, being a total fluence ~3X10% erg cm? (15-150 


keV), the second highest value ever recorded by Swift!°, and c) by the faintness of the X-ray 
counterpart, as the X-ray flux to gamma-ray fluence ratio, log fx,11nr/Fy ~ -7.9, lies below the typical 
GRB distribution! as expected for an explosion in a rarefied circumburst medium!” (Methods). 
The association with a galaxy at 346 Mpc implies that GRB211211A is one of the closest long 
bursts ever discovered, yet the properties of its gamma-ray emission, such as the negligible 
temporal lag, short variability timescale, and hard spectrum, do not fit into this class of events 
(Extended Data Figure 2). These are distinctive features of short bursts and classify GRB211211A 
as a hybrid event, analogous to GRB060614°. In addition to its prompt gamma-ray phase, several 
lines of evidence differentiate GRB211211A from canonical long GRBs. The GRB does not lie in 
a star-forming region (Methods) and late-time optical imaging rules out any bright supernova (SN) 
at its location (Extended Data Figure 4): as the dust content along the line of sight is negligible, a 
luminous SN similar to SN 1998bw'* is excluded out to z~0.8. A faint and short-lived SN similar 
to SN2008ha° is also ruled out by the optical limits. The GRB location and the global properties 
of its host galaxy provide indirect evidence for a stellar progenitor different from a collapsing 
massive star and are instead consistent with a binary NS merger (Methods). 

The unambiguous proof of a compact object binary merger comes either from its gravitational 


wave signal? or from its kilonova, a short-lived glow of quasi-thermal radiation powered by the 


radioactive decay energy of heavy nuclei!*, produced in the merger ejecta via rapid neutron capture 
process (r-process). The first known kilonova was AT2017gfo, characterized by an early (<12 
hrs) UV/optical peak'® followed by a longer-lasting infrared signal'’!*!°. We find that a similar 
component is identified in the UV/optical/infrared (UVOIR) counterpart of GRB211211A, 
providing us with the direct link to compact binary mergers. 

The multi-wavelength emission that follows a GRB is the superposition of multiple components. 
The dominant one is usually the afterglow, a broadband synchrotron radiation emitted by a 
population of electrons shock-accelerated by the GRB outflow7’. We use the X-ray data to probe 
the contribution of this non-thermal component. The X-ray spectrum is well-described by a power- 
law with slope 6x~0.5 and negligible absorption along the line of sight. When extrapolated to 
lower energies, this model roughly matches the observed optical fluxes at To+1 hr and shows no 
evidence for an additional component at this time. However, at later times, the multi-frequency 
spectral energy distribution (Figure 2) identifies emission in excess of the standard afterglow: the 
UVOIR counterpart is consistently brighter than the extrapolation of the non-thermal power-law, 
and is characterized by a steeper spectral index Buvoir > 2 for t> 1d. Its spectral peak lies in the 
UV range (u band, observer-frame) at To+0.2 d and then progressively cools down to near-infrared 
(nIR) wavelengths (K band at ~To+4 d). 

We rule out that a reverse-shock (RS) powered afterglow or a supernova (SN) onset could explain 
this low-energy component. The former arises within the GRB outflow and is characterized by an 
optical rebrightening peaking from a few seconds to ~1 hr after the burst”!??. However, RS 
emission quickly cools off and shifts to the radio band, typically within the first day after the burst. 
This is not consistent with the observed SED evolution. Moreover, a low nickel-yield explosion 


would also produce a short-lived UV/optical flare powered by shock heating in the SN blast- 


wave’>. We studied a broad range of collapsar-associated SNe varying nickel-yields, stellar 
properties, and explosion energies. Although this model explains the lack of a bright SN at late 
times and can reproduce the basic features of the early optical emission such as the bolometric 
luminosity and photospheric radii, the predicted spectrum is too hard (Extended Data Figure 6): 
UV emission dominates and we cannot reproduce the bright and long-lived nIR emission without 
the addition of a second, neutron rich outflow (see Methods). 

After subtracting the afterglow contribution from the data, we find that the UVOIR excess is well 
described by a thermal spectrum and that the best-fit parameters point to a hot (7~15,000 K, rest- 
frame) fireball in rapid expansion with apparent velocity v=0.5 c. The luminosity, temperature and 
emitting radius of this thermal component display a striking resemblance to AT2017gfo~4 (Figure 
2), and we interpret it as the kilonova emission associated with GRB211211A. A kilonova in 
GRB211211A and consequently its association with a compact binary merger tie the lack of SN, 
the GRB environment, and the evolution of its UVOIR counterpart in a coherent explanation. 
Our dataset allows us to probe the earliest phases of the kilonova onset, not observed in the case 
of AT2017gfo. Although the broadband emission is initially dominated by the non-thermal 
afterglow, evidence for a thermal component is found as early as To + 5 hr. Figure 3 shows the 
different behaviors of the X-ray and UVOIR counterparts. The latter requires an additional 
component, which we model using simulated kilonova light curves”> with wind ejecta mass Mw in 
the range 0.01 - 0.1 Msun, and dynamical ejecta mass Mayn~0.01-0.03 Msun. The ejecta velocity and 


kilonova bolometric luminosity, Lpor~ 3X10” erg s"! (isotropic equivalent) inferred at early times, 
y g pic eq y 


are challenging to reproduce with purely radioactive-powered models”, even when accounting for 
different density profiles and the larger projected area along the polar axis*® (Methods). We 


therefore explore alternative models in which the merger ejecta is re-energized by a central engine 


or modified by the interaction with the GRB jet. The former group of models, envisioning either 
a highly magnetized NS or fallback accretion onto the central black hole (BH), is often invoked to 
explain a long-lasting gamma-ray emission”’?*. However, an active engine would leave observable 


t??, which are not consistent with its timescales (too early) or colors 


imprints on the kilonova ligh 
(too red) (Extended Data Figure 6). 

We therefore consider a model in which jet-ejecta interactions shape the observed emission. As 
the jet propagates through the massive (>0.01 Msun) cloud of radioactive ejecta, it heats and 
partially disrupts its density structure, carving a funnel of low-opacity low-density material along 


the polar axis*’. By exposing the inner, hotter surface of the ejecta, an energetic (E),iso~6X 10°! erg) 


GRB jet makes the kilonova emission both bluer and brighter*° for an observer close to its axis. 
Shock-heating may also contribute to distribute the energy. Viceversa, the ejecta imparts a wide 
angular structure on the GRB jet before it breaks out?!3*. High latitude emission from the jet wings 
arrives at a later time because of the longer path that photons travel and may produce a low- 
luminosity fast-fading X-ray transient** consistent with the observed X-ray behavior. This feature 
may become visible in the case of a “naked” structured GRB jet expanding into a low-density 
circumburst medium. 

We conclude that, although the long duration of the prompt phase challenges our understanding of 
compact binary merger models, a merger progenitor naturally explains all the other observed 
features of GRB211211A. At a distance of 346 Mpc, this GRB lies within the distance horizon of 
forthcoming gravitational wave (GW) observing runs*4 and, had the GW network been online at 
the time of the burst, this event would have likely resulted in a joint detection of GWs and 
electromagnetic (EM) radiation. We note that its EM properties are very different from the multi- 


messenger transient GW170817: whereas the EM counterparts of GW1708177'>!*!7 would be 


challenging to localize beyond ~150 Mpc, GRB211211A would be visible out to z~1 by most 
space-borne gamma-ray detectors. Moreover, rapid X-ray and UV/optical follow-up of GW 
sources would detect its counterpart out to z~0.2 assuming a sensitivity comparable to Swift. 

To determine the rate of hybrid GRBs, we examine the Swift GRB catalog'® in search of bursts 
similar to GRB211211A and GRB060614. At large distances (z~1), their classification would rely 
solely on the high-energy properties which point to regular bursts of long duration (Extended Data 
Figure 7). Without a systematic study of GRB lags, spectra and durations, it is not possible to 
assess the total number of hybrid bursts thus far detected. Therefore, we turn to lower redshifts 
where a clear observational signature of these events is the lack of a SN. SNe associated with 
GRBs! peak between My~-18.5 mag and My~-20 mag, and sensitive SN searches are regularly 
undertaken for GRBs within a redshift z<0.3, which we identify as the maximum distance for a 
homogeneous identification. We review the entire sample of Swift bursts with duration >2 s and a 
putative host galaxy at z<0.3 and find a total of 20 GRBs in 17 years of mission (2005-2021). Of 
these, 8 are associated with a SN, 3 have no meaningful constraints, and 9 have deep limits on any 
accompanying SN. The chance alignment between a bright galaxy and an afterglow with sub- 
arcsecond localization is typically <1% (ref. 11), thus it is unlikely that all nine bursts are distant 
background objects. Furthermore, four of them (GRB060614, GRB060505, GRB191019A and 
GRB211211A) have UV counterparts constraining their distance scale>*>3°. We conclude that at 
least some of these long duration bursts are physically associated with a low-redshift galaxy and 
lack a SN, forming a new class of hybrid GRBs produced by compact binary mergers. After 
accounting for Swift’s duty cycle and limited field of view, we derive a volumetric all-sky rate of 
0.04-0.7 Gpe? yr! (68% c. 1.; ref. 37), lower than the local rate of short GRBs**. The true rate of 


events depends on the unknown beaming factor of these outflows, but, assuming similar jet 


properties to short GRBs*’, hybrid long duration bursts may represent ~10% of the population of 


EM counterparts to GW sources caused by compact binary mergers. 
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Figure 1 - The field of GRB211211A 

This false color image (a) combines optical (blue) and near-infrared (red+green) HST observations 
of GRB211211A. Two bright galaxies (G1 at z~0.076, and G2 at z~0.458) and several fainter ones 
are visible, but no source is detected at the location of GRB211211A. The most likely host galaxy 
is G1, a low mass late-type galaxy. The projected physical offset between the burst and the galaxy’s 
center is ~8 kpc, one of the largest ever measured for a long burst. The same field is shown in the 
UV w2 filter (b) and optical I filter (c) at 1 hr and 10 hr after the burst, respectively. The solid lines 
show the slit position used for optical spectroscopy with Gemini/GMOS-S. The bright UV 
counterpart rules out a high-redshift origin, whereas its rapid reddening is consistent with the onset 


of a kilonova. 
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Figure 2 - Spectral evolution of the GRB afterglow and kilonova 

The spectral energy distribution (a) combines gamma-ray (diamonds), X-rays (circles) and UVOIR 
(squares) data at different times, as indicated by the labels. It shows that non-thermal radiation 
(solid line) dominates at early times and at higher energies. At lower energies, we identify the 
emergence of a thermal component peaking at blue wavelengths at 5 hr, and rapidly shifting toward 
redder colors. Downward triangles are upper limits (3 oc.1.). For plotting purposes, each epoch 
was rescaled by the following factors (from top to bottom): 1, 1, 10°°8, 10°!°, 10°74, 10°. The 
bolometric luminosity (b), temperature (c), and emitting radius (d) of the thermal component are 
similar to AT2017gfo*4 (gray circles). Solid lines show the best fit power-law models to the 
dataset. Dashed lines in panel d show the predicted radius for constant expansion velocities of 0.3 
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Figure 3 -A kilonova in the long GRB211211A 

Multi-color light curves in X-rays, UV (uvm2), optical (BRI) and infrared (K) are compared to 
models predictions of a kilonova (solid line) in addition to the non-thermal emission (dashed line). 
The shaded area shows the range of possible fluxes reproduced by kilonova simulations with wind 
mass My between 0.01 Mgun (lower bound) and 0.1 Msun (upper bound), and dynamical ejecta mass 


Mayn between 0.01 Msun (lower bound) and 0.03 Msun (upper bound). 
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Extended Data Figure 1 - Prompt gamma-ray phase of GRB211211A 

The Swift background subtracted light curves of GRB211211A are shown in two energy bands, 
and compared with the time history of GRB060614 (gray shaded area) rescaled at a distance of 
346 Mpc. The time bin is | s. Both bursts display a first episode with hard spectrum (dominant in 
GRB211211A), followed by a long-lasting tail with soft spectrum (dominant in GRB060614). The 
inset zooms in the first 12 s, showing a weak precursor at To preceding the main prompt event. The 


time bin is 16 ms. 
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Extended Data Figure 2 - GRB classification scheme 

The traditional GRB classification, based on the duration/hardness ratio diagram (a), is not 
unambiguous. Additional classificators, used to break the degeneracy, are the lag-luminosity 
relation (b), the variability timescale (c), and the Amati relation (d). Long GRBs (circles) and short 
GRBs (squares) occupy different regions of these plots. Dashed lines show the boundaries of the 
long GRB regions (shaded areas). GRB211211A (star symbol) belongs to the class of long soft 


bursts (a), but its other high-energy properties are common among short GRBs. 
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Extended Data Figure 3 - Constraints on the distance scale of GRB211211A 

The nIR brightness of GRB host galaxies (short GRBs: squares; long GRBs: circles) is reported as 
a function of their redshift. For comparison, a randomly selected sample of field galaxies from the 
CANDELS survey is shown in the background (octagons, with the symbol size proportional to the 
galaxy mass). The non-detection of an underlying galaxy in deep HST F/60W imaging rules out 
most of the parameter space occupied by GRB hosts (hatched area). Additional constraints from 
the UV and X-ray afterglow rule out the case of a GRB in a distant (z>1.5; shaded areas) faint 
galaxy. These observations support the physical association between GRB211211A and the nearby 


galaxy at z=0.0762 (star symbol). 
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Extended Data Figure 4 - No supernova associated with GRB211211A 

Optical upper limits in the r-band (red) and i-band (blue) rule out the presence of any known 
supernova following GRB211211A. Bright SNe associated with GRBs, such as SN1998bw!4 and 
SN2006aj**, would have been detected up to z=0.8 (dashed line) and z=0.65 (dotted line), 
respectively. Symbols show the peak magnitude of core-collapse SNe from the ZTF Bright 
Transient Survey Sample®’ rescaled at z=0.3, demonstrating that most ordinary SNe were 
detectable up to this distance. At the distance z~0.076 of the putative host galaxy, the faint 


SN2008ha‘ (solid line) is also ruled out. 
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Extended Data Figure 5 - Host galaxy spectral energy distribution 


(a) The model SED (blue line) and model photometry (blue squares) derived using Prospector is 
compared to the observed photometry (red circles). The fit residuals (b) are displayed in the bottom 


panel. The inset (c) displays a Gemini/GMOS-S spectrum of the host galaxy in the vicinity of Ha, 


[NII], and the [SI] doublet, yielding z~0.0762. 
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Extended Data Figure 6 - A comparison of models for luminous blue transients 
The UVOIR counterpart is compared with a set of models producing luminous (Lpor~10* erg s*!) 
and short-lived (< 7 d) transients: (a) a low-nickel SN from a fallback collapsar?? underpredicts 
the optical/nIR emission; (b) a magnetar-powered kilonova”® does not reproduce the timescales 


and colors. 
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Extended Data Figure 7 - GRB211211A at high redshift 
Simulated Swift light curve in the 15-150 keV energy band (observer’s frame) of GRB211211A 
assuming z=1, an incident angle of 45 deg, and ~6,000 cts s! background rate. Time bin is 1 s, 
error bars are 68% confidence level. An event similar to GRB211211A would be detected up to 
redshift z~1 and beyond. At these distances, it would appear as a standard long GRB with a 


duration of ~20 s. 


Methods 


Observations and Data Reduction 

Swift - The Swift X-ray Telescope (XRT) began observations of GRB211211A at 69 s after the 
GRB trigger and continued monitoring for ~3 days. The X-ray lightcurve and spectra were 
retrieved from the Swift XRT GRB repository*’. The GRB counterpart was detected by the 
UltraViolet and Optical Telescope (UVOT) in all filters. UVOT data were analyzed using 
HEASOFT v.6.30. Individual exposures from the same epoch were aligned and then combined 
using the task uvotimsum. To minimize the contamination from the nearby galaxy Gl, we 
determined the source count-rates in 3” circular apertures with the background estimated in a 
concentric source-free annular region. We used an aperture correction to determine the equivalent 
source count-rates in a 5” aperture, and then applied standard UVOT magnitude zeropoints*!. Late- 
time images were used to determine the host galaxy brightness using a 5” circular aperture. The 
photometry of the transient and host galaxy are reported in Extended Data Table 1. 

Optical/nIR imaging - Our observations of GRB211211A began at To+5.1 hrs with the 50-cm 
Multicolor Imaging Telescopes for Survey and Monstrous Explosions (MITSUME*34) 
acquiring 11060 s simultaneous images in g, R., and J. filters. Observations were carried out for 
three consecutive nights, but the GRB counterpart was only detected on the first night. Data were 
reduced using the MITSUME pipeline* and image subtraction between the first, second and third 
epoch was performed using HOTPANTS**. Upper limits were derived by planting artificial 
sources of known brightness, then repeating the image subtraction step. 

Starting on 2021 Dec 11 UT 22:11:29 (~To + 9.0 hr), multi-band (U, B, V, R, J) optical observations 
of GRB211211A were carried out using the 4Kx4K CCD Imager mounted at the axial port of the 


3.6m Devasthal Optical Telescope (DOT*’**4?>°). Observations continued until twilight time 


(~To + 11.4 hr) and were repeated in subsequent nights until the optical counterpart was no longer 
detected. Starting on 12 Dec 2021 at 3:53 UT (To + 15 hr), additional imaging of the GRB field 
was carried out with the 2.2m telescope, equipped with the Calar Alto Faint Object Spectrograph 
(CAFOS), at the Calar Alto Observatory in Almeria (Spain) using the gr’ Sloan filters. 
Additional epochs of observations were performed between 2021 Dec 31 and 2022 Jan 02 to search 
for the possible SN peak. We acquired multiple R-band images (24X300s) with the 1.3m Devasthal 
Fast Optical Telescope (DFOT) and deep g'r'i’ exposures with the CAHA telescope. 

At late times, we targeted the field with the Large Monolithic Imager (LMI) on the 4.3-m Lowell 
Discovery Telescope (LDT) using the u, i, and z filters in order to refine the photometry of its host 
galaxy. 

The data were pre-processed using standard CCD reduction techniques including bias subtraction, 
flat-fielding, fringe correction and cosmic ray rejection. Aperture photometry on the GRB 
counterpart was performed using Source Extractor?! and selecting circular apertures with radius 
1.5 times larger than the image full width half maximum. Forced photometry was performed on 
images acquired at later times (> To + 2 d). Host galaxy magnitudes were derived using Kron-like 
elliptical apertures (MAG AUTO). The photometric zeropoints were calibrated using the same 
nearby point sources in the Sloan Digital Sky Survey (SDSS*’) and transformed using empirical 
equations*’. 

We imaged the field with the Hubble Space Telescope (HST) using the WFC3/UVIS and IR 
cameras with the FS$/4W and F/60W filters, respectively. The data were reduced using standard 
procedures within the DrizzlePac software™ in order to align, drizzle, and combine exposures. We 
used Source Extractor to detect sources and perform aperture photometry. For each filter, we used 


the zeropoints stored in the keyword PHOTFLAM. 


Gemini Spectroscopy - We utilized the Gemini Multi-Object Spectrograph (GMOS) mounted on 
the 8.1-m Gemini South telescope to obtain a series of 4X600 s spectra, using the R400 grating 
and a 1” slit. We chose a slit position angle (PA) of 64 deg such that the slit covers both the host 
galaxy’s center and the GRB’s optical position. The data were reduced and analyzed using standard 
procedures in Gemini IRAF. In order to correct for slit losses, the flux-calibrated spectra were 
matched to the photometry of the host galaxy. Line fluxes were derived by fitting each line with a 


Gaussian function and estimating the continuum from nearby spectral regions. 


Classification of GRB211211A 

GRBs are classified based on the properties of their prompt gamma-ray phase. The prompt 
emission of GRB211211A (Extended Data Figure 1) displays three different episodes: a weak 
precursor, a bright multi-peaked main burst, and a highly variable temporally extended emission. 
The time intervals for spectral and temporal analysis were selected to characterize them separately. 
Swift and Fermi data were processed using HEASOFT v.6.30. Spectra were extracted from the 
Fermi Gamma-ray Burst Monitor data and fit within XSPEC. The temporal properties were 
derived from the Swift BAT light curves using well-established techniques**>’. 

The precursor phase has a short duration of 0.15 s, a soft spectrum peaking at ~75 keV, a minimum 
variability timescale of 21+4 ms, and a significant lag 13:= 1613 ms between the temporal 
structures observed in the 50-100 keV (band 3) and in the 15-25 keV (band 1) energy bands, 
respectively. At 346 Mpc, the measured flux of 8x10°’ erg cms! (10 - 1,000 keV) corresponds 
to a luminosity of only ~10* erg s‘!. 


After a 1 s period of quiescence, we detect the onset of the main prompt emission, which consists 


of multiple overlapping peaks lasting for approximately 10 s. The time-averaged spectrum peaks 


at 750+10 keV, the minimum variability timescale is 14+5 ms, and the temporal lag is negligible 
with t3:= —0.9*23 ms. The total fluence measured during this episode is ~3.7x 107 erg cm? (10 - 
1,000 keV), one of the highest ever measured for a GRB. However, for a distance of only 346 Mpc 
the total isotropic-equivalent gamma-ray energy Eyiso would be ~5X10°! erg within the typical 
GRB range”. 

A brief (3 s) period of low-level persistent emission precedes the onset of a long-lasting tail. The 
time-averaged spectrum of the extended emission has a soft peak of 52+2 keV, the minimum 
variability timescale is 42 +9 ms, and the lag, t31= 7*3 ms, is positive. The total fluence is ~5x 10° 
erg cm (10 - 1,000 keV), corresponding to Eyiso ~7X10°° erg. 

We compare these properties to the population of GRBs using four classifiers: the 
duration/hardness-ratio diagram>”, the lag-luminosity relation”®, the variability timescale>’, and the 
Amati correlation®® (Extended Data Figure 2). Similarly to GRB060614°, GRB211211A shows 
characteristics that are intermediate between the two main GRB classes: the traditional 
classification based on duration and hardness ratio places this event in the class of long GRBs, 
however its other properties fit within the class of short bursts. Its hybrid nature does not allow us 


to unambiguously link it to a progenitor system based solely on its high-energy properties. 


The GRB environment and its host galaxy 

The GRB environment typically offers stringent, albeit indirect, evidence of its progenitor system. 
In the case of GRB211211A, no underlying host galaxy is detected in late-time HST imaging 
(Figure 1). By planting artificial sources with an exponential disk profile and different brightness, 
we derive an upper limit of F$/4W>26.5 AB mag and F'160W>27.6 AB mag. Since no coincident 


galaxy is found, we analyze the GRB field to search for its most likely host. We identify 7 galaxies 


within 10’’ from the GRB position (see Figure 1): G1 with r = 19.50 +/- 0.02 mag at an offset of 
5.55”’ +/-0.03’’, G2 with r = 20.88 +/-0.05 mag at an offset of ~10’’, and 5 faint (x > 26 AB mag) 
extended objects at an offset between 2.5’’ and 10°’. By using the galaxy's number counts in r- 
band®°, we derive a chance alignment of 1.4% for G1, >10% for G2, and >40% for the other faint 
galaxies. Therefore, probabilistic arguments favor the association between GRB211211A and G1. 
Moreover, in our spectroscopic observations we find no evidence for any emission lines at the 
GRB position down to > 2 x10!’ erg cm? s"! A"! in the range 4,800-6,100 A. Using [OII] 3727 
and Hf as indicators of unobscured star formation®', we place an upper limit on the star formation 
rate SFR<1 Mg yr’ for z<0.65. This corresponds to the median SFR of long GRB hosts” at z < 1, 
providing additional constraints on any possible underlying galaxy. 

The spectrum of G1 shows several emission lines including Ha, [NII], and [SII] at a common 
redshift of z=0.0762 +/- 0.0003, consistent with the previous report? based on data from the Nordic 
Optical Telescope (NOT). Assuming a ACDM cosmology® with Ho=69.8 km Mpc"! s, we find 
a luminosity distance di = 346 Mpc, and a distance modulus y= -37.7 mag. Using the host galaxy 
photometry (Extended Data Table 1), we compute a rest-frame absolute B-band magnitude of 
Mpx-17.6 AB mag, corresponding to Lg~0.1L* when compared to the galaxy luminosity 
function™ at a similar redshift (0.05<z<0.2). 

The brightness (Lya~10*° erg s‘') and relative ratio of these lines (log [NII]/Ha ~ -0.7) point to a 
star-forming galaxy with SFR~0.05 Mo yr’ and sub-solar metallicity 12 + log O/H ~ 8.4. We also 
find evidence for weak [Mg I 45175A] absorption at ~5567 A, indicative of an evolved stellar 
population, although this feature is affected by a nearby sky line. 

We model the galaxy’s surface brightness using GALFIT®. A good description (72~1.03) of its 


morphology is obtained by including two Sersic profiles with index n=1, one with half-light radius 


Re1~2.15 arcsec (F814W; ~3.1 kpc at z=0.076) and another one with Re,~0.5 arcsec (F814W; ~0.7 
kpc at z=0.076) to model the central bar. Similar results are obtained on the F/60W image with 
Rei~2.34 arcsec and Re ~0.64 arcsec. The half-light radius rso ~ 1.1 arcsec obtained through 
Source Extractor is given by the weighted average of these two components. 

The galaxy’s global properties were determined by modeling its SED (Extended Data Table 1) 
with Prospector, adopting the same settings used for GRB host galaxies!!°’. We derived a stellar 
mass of M= 0.9734 x 10° Mo, a star-formation rate SFR=(0.06 + 0.02) Mo yr, a low dust content 
Ay = 0.0910:38 mag, and a mass-weighted stellar age t = 5*3 Gyr. When compared to the sample 
of long GRBs, the properties of the host of GRB211211A are not unprecedented but extremely 
uncommon. The inferred SFR lies in the bottom 10% of the observed distribution, leading to an 
unusually low specific SFR, sSFR~0.06 Gyr!. This value is below the main sequence of star 
forming galaxies®, indicating that G1 may be migrating to a quiescent phase. This differs from 
the typical environment of long GRBs at both high and low-redshifts: for comparison, nearby 
events such as GRB060218 and GRB100316D were associated with sSFR~4 Gyr! and 
sSFR~0.2 Gyr!, respectively’. Dissimilarities with the class of short GRBs also exist: the 
stellar mass lies at the bottom 10% of both short GRB and SN Ia host galaxies’!’”, as for 


GRB060614 which was hosted by a dwarf galaxy°. 


Spectral Energy Distribution 
The spectral energy distribution (SED) of the GRB counterpart at different times is shown in 
Figure 2. These epochs were selected to maximize simultaneous multi-wavelength coverage. 


When needed, the data were rescaled to a common epoch using the best fit temporal model. 


In the first epoch at To+100 s, the X-ray emission is characterized by a flat spectral index Bx= 0.00 
+/- 0.03. A spectral break is required above ~10 keV to account for the lower flux and soft spectral 
index, Bpat~2, measured in the hard X-ray band. In addition, the high X-ray-to-optical flux ratio, 
Fx/Fo*100, requires a turn-over to a steep spectrum between the X-ray and optical band. These 
properties are consistent with self-absorbed synchrotron radiation in the fast cooling regime. The 
location of a self-absorption frequency, va~10 eV, indicates a compact emitting region’? with 
radius R~10'? ([/300)*4 cm, where I is the outflow bulk Lorentz factor. This radius is typical of 
dissipation processes within the GRB outflow, indicating that at ~To+100 s the prompt phase is 
still dominant at both X-ray and optical wavelengths. 

In the second epoch at To+1 hr, the GRB counterpart displays blue colors with a spectral index Bo 
= 0.23+0.10 in the UV and optical bands. At X-ray energies the spectrum, extracted between 3 ks 
and 5 ks, has a slope of Bx = 0.50+0.05. This index points to synchrotron radiation in the slow 
cooling regime, in which the cooling frequency is ve> 10 keV and the synchrotron frequency is 
VmS1 eV. In this case, the X-ray spectral slope is related to the energy distribution of the emitting 
electrons, N(E)« E? with p = 2 Bx + 1 = 2.0+0.1. This is a fundamental constraint to the long- 
term afterglow evolution. The steepest spectral slope explained by this model is p/2~1.05, and 
only for energies above vc. Therefore, the UVOIR and X-ray non-thermal afterglows are bound to 
remain on the same spectral segment for times >To+1 hr. 

Starting from ~To+5 hr, a simple non-thermal spectrum can no longer reproduce the broadband 
emission. An UVOIR excess is detected at all epochs. It is characterized by a narrow spectral shape 
peaking in the uw band, well described by a blackbody function with temperature 7~15,000 K (rest- 
frame) and a luminosity Lyo® (3.542.0) x 10” erg s‘!. We therefore fit each SED epoch with a 


blackbody (UVOIR) plus power-law (X-ray) model, and derive the total integrated blackbody 


luminosity, its temperature and radius as a function of time (Figure 2). The luminosity is better 
constrained in our second epoch at To+ 10 hr, Lio: = (1.940.2) x 10% erg s"! and is seen to decrease 
in time following a power-law « f°>, consistent with the evolution of AT2017gfo**. 

GRB distance scale 

We investigate the joint X-ray/UV/optical SED at 1 hr to place a direct upper limit on the GRB 
distance scale. UVOT spectra were created with the tool wvot2pha using the same source and 
background regions selected for photometry. We adopt a power-law model and include the effects 
of absorption, dust reddening, and intergalactic medium attenuation as implemented in the XSPEC 
models zphabs, zdust, and zigm. The Galactic absorption was fixed to Nu = 1.76 x 107° cm”? and 
the reddening at E(B-V) = 0.015 mag. All other parameters were left free to vary. We increase the 
redshift from 0 to 2.5 in steps of 0.1 and find the best fit model by minimizing the Cash statistics, 
recording its value at each step. Based on the variations of the test statistics, we derive an upper 
limit of z < 2.3 (99.9% confidence level) from the UV/optical data, and z< 1.5 (99.9% confidence 
level) from the joint X-ray/UV/optical fit. By imposing the redshift of the putative host galaxy, 
z~0.0762, we find no evidence for any dust extinction or absorption at the GRB site with 3 o upper 
limits of E(B-V). < 0.005 mag and Nu, <9 x 10!% cm”, respectively. This is consistent with the 
location of the GRB, well outside the galaxy’s light. 

Origin of the X-ray afterglow 

Swift observations show a rapidly fading X-ray afterglow followed by a shallower decline Fx r@ 
with a=1.1149-38 between | ks and 40 ks, and a final steep decay with a=3 + 0.5 after 40 ks. 
Based on this model, we infer an X-ray flux of ~4 x 10° erg cm” s"! at 11 hrs. This corresponds 
to a luminosity Lx~6x10* erg s"! at a distance of 346 Mpc , nearly two orders of magnitude below 


the typical X-ray luminosity of cosmological GRB afterglows at this epoch (cf. Figure 7 of ref. 21). 


The low ratio between the observed X-ray flux and the emitted gamma-ray fluence, log fx,11nr/Fy 
~ -7.9, is indicative of atypical properties for this explosion (cf. Figure 17 of ref. 11). 

Our SED analysis (Figure 2) demonstrates that the X-ray counterpart is dominated by non-thermal 
emission consistent with synchrotron radiation. While we interpret the early (<300 s) X-ray 
emission as the tail of the prompt phase, at later times (>1,000 s) the most common origin of non- 
thermal afterglow radiation is the interaction between the ambient medium and the GRB jet 
occurring at large distances (>10!’ cm) from the central source. In this external shock model”, a 
flux decay rate of 2 or faster is explained by geometrical factors due to the collimation of the GRB 
outflow’. The time t;at which the light curve steepens, the so-called jet-break, increases with the 
jet opening angle @.. A jet-break at 40 ks would require a very narrow Jet, and then can only achieve 
a decay of a=p~2.1, in mild tension with the observations. We tested the hypothesis of an early 
jet-break by modeling the X-ray and early (~To+1 hr) UVOT data with afierglowpy”® assuming a 
uniform external environment and both a top-hat and a Gaussian lateral structure for the jet. 
Despite the dataset being limited, it provides tight constraints to the model: the flat UVOT SED at 
To+1 hr (Figure 2) requires the synchrotron peak to lie close to the optical range, constraining the 
value of vm and F’px; the X-ray spectrum places the cooling frequency at v-> 10 keV and provides 
a measurement of p~2.0-2.1, and the X-ray light curve constrains the jet opening angle 0. and the 
viewing angle 6,. We performed Bayesian parameter estimation with emcee” and nine free 
parameters: Ex, iso, @, Oy, an outer jet truncation angle 6, 1, p, &e, €g, and the participation fraction 
Ex. The best fit has a reduced chi-squared y?~1.8, fits with En frozen at 1 found similar y? but 
required unphysical shock parameters ¢- © eg © 1. The parameter estimation reports a jet of energy 
Ek,iso ~ (0.8-17) X10°! erg, width @ ~ 1.9-5.7 deg, viewed at 0, ~ 1.1-5.4 deg from the jet axis. 


The external density is n ~ 0.016 - 12 cm. The shock parameters are p ~ 2.1-2.2, & ~ 0.05 - 0.77, 


ép ~ (0.1 - 6.0) X10%, and En ~ (0.8-9.6) 10°. The beaming-corrected kinetic energy of the jet in 
this scenario is 0.4 - 4.4 x10? erg. Assuming the angular size-corrections between the afterglow 
and prompt emissions are similar, this scenario gives ~65% probability to an unphysical gamma- 
ray efficiency ny = Ey,iso/Ex,iso > 100% and a 90% probability ny > 15%. In combination with the 
poor reduced chi-squared of 1.8 we conclude it is challenging for an external shock to 
simultaneously reproduce the salient features of the GRB afterglow: a flat UV/optical spectrum at 
Tot+l hr, an X-ray spectrum (x ~ 0.5, and a steep decay of the X-ray flux after 40 ks, while 
remaining within the energetic limits of the prompt emission. 

If not produced by a jet break, a rapid drop in brightness is difficult to produce due to the relativistic 
and extended nature of the GRB outflow. Due to the curvature effect!?-’* any rapid decrease in 
brightness in the lab-frame of the GRB will be smeared out in the observer frame due to different 
arrival times of the photons, producing a decay of a=2+Bx~2.5. Nevertheless, this is a steeper slope 
than that allowed by the jet-break model and may present a better description than the standard 
external shock. If interpreted as a curvature effect, the steepening at 0.5 d links the observed X- 
ray emission either to long-lasting activity of the central engine, as in the “internal plateau” 
model’?®°, or to the angular structure of the GRB jet. If a structured jet produces GRB prompt 
emission in the high latitude regions (the jet “wings”’), this emission would be Lorentz de-boosted 


relative to the core prompt emission and delayed via the curvature effect** 


. With appropriate jet 
structures, this can manifest as X-ray emission with a shallow decay followed by a steep declining 
light curve. This feature, normally hidden by the brighter external shock emission, may become 
apparent in the case of a “naked’’ structured GRB exploding in a rarefied medium. This latter 


model offers a consistent explanation of the X-ray behavior of GRB211211A and its physical 


offset from the galaxy without the requirement of hours-long activity of the central engine. 


Despite uncertainty in the physical origin of the afterglow emission, the observed X-ray spectrum 
is well measured and its extrapolation to the UVOIR bands unambiguously places it below the 
UV/optical detections after ~To+5 hr. The observed UVOIR excess was measured by subtracting 
this extrapolated non-thermal component. This procedure does not require a physical 
interpretation of the non-thermal emission and provides an upper bound on the non-thermal 
contribution in the UVOIR bands. Thus, the identification of the UVOIR excess does not depend 


on the specific physical interpretation of GRB211211A’s non-thermal emission. 


Origin of the UVOIR excess 

Collapsar model - We first examine the most common case of a long GRB produced by the 
collapse of a rapidly rotating massive star (collapsar). The emergence of the SN blast wave can 
produce a luminous blue emission in excess of the standard afterglow”?, and we test whether this 
is consistent with the observed UVOIR excess in GRB211211A. Collapsars arise from compact 
stellar cores and produce energetic and long-lived Ic supernovae/hypernovae. However, if the 
collapsar engine does not produce significant *°Ni (e.g. from a fallback collapsar), the blastwave 
produces a short-lived SN light curve that dies out in the first 10 d. To test this model, we ran a 
series of hypernova explosions, varying the mass (2.5-40 Mj) and density profile (varying the 
slope in the density of the core and envelope) of the progenitor star as well as the explosion energy 
(spherically 10°!-10° erg). Although we can reproduce the evolution of the bolometric luminosity, 
the early-time emission in our best-fit model is too energetic (in the UV and extreme UV). As the 
ejecta cools, the emission peaks in the IR at late times, but the luminosity is several orders of 
magnitude too dim to explain the observations. To account for the optical and IR emission, the 


photosphere of the rapidly expanding SN must uncover the collapsar accretion disk and wind ejecta 


from this disk must have similar-enough properties to NS merger disks*!** to produce a kilonova- 
like transient. For the collapsar model to work, we must also explain the offset of the GRB from 
its host galaxy. O/B stars in binaries can be unbound during the supernova explosion of the 
primary star, imparting a “kick” of up to 200 km s"! onto the O/B companion®. This proper motion 
could move the companion O star well beyond its star forming region (~1 kpc in 5 Myr), but it is 
unlikely that this kick is sufficient to explain the large offset of this burst. In summary, a massive 
star progenitor for GRB211211A would naturally account for its long duration but requires a 
combination of unusual circumstances (a low °°Ni yield explosion, a neutron rich disk outflow, 


and an extreme kick velocity) to explain the entire set of observations. 


NS merger model - The observed excess emission is much better fit by the disk-wind ejecta from 
a NS merger. Figure 3 shows the range of model predictions consistent with the observations: only 
a small subset of light curves (4 out of 900 in the “on-axis” angular bin; 6,~0-16 deg) match the 
observing constraints. The bright UV/optical counterpart suggests a massive (>0.01 Ma) wind 
component to the kilonova ejecta. However, the time-dependent spectra from the Los Alamos 
National Laboratory (LANL) grid of kilonova models*> produce light curves that are too dim to 
match the observed UV/optical luminosities or require too large an ejecta mass (~0.1 Mc). Models 
with large ejecta mass (M,,~0.1 Mo) fit better the early time data but overpredict the fluxes at later 
times and, vice versa, the model with lower ejecta mass (M,~0.01 Me) provides a good description 
of the dataset only after ~11 hrs. All consistent models adopt a toroidal morphology for the high- 
opacity ejecta and a polar outflow of low-opacity ejecta and high expansion velocity vw~0.3 c. 

It is likely that a number of alterations to the kilonova ejecta mechanism can help explain the early 


excess emission. For example, we have not conducted a detailed study varying the composition 


that changes both the opacity and the radioactive heating. Uncertainties in radioactive energy 
deposition* and in the properties of the disk-wind ejecta allow for a wide range of behaviors and 
our study here only touches the surface of all possibilities. However, in its simplest form, a 
radioactive-powered kilonova captures the late-time evolution of the observed UVOIR transient 
but has difficulties in reproducing the bright optical emission seen at early times (To+0.2 d). 

An alternative way to alleviate the requirement on the ejecta mass is that the kilonova is powered 
by an additional energy source or affected by the jet-ejecta interactions*®. To study the engine- 
powered models, we used the same method as in previous studies*®. For central power sources, 
either a magnetar or fallback accretion on the central BH, the energy must transport out from the 
center to affect the light-curves. In these models”*, the central power sources do not alter the 
emission until ~5 d after the merger for wind mass ~0.01 Mo. However, if the jet is able to evacuate 
a region above the compact remnant, this delay can be reduced. We mimicked this evacuation by 
a series of spherically symmetric models, reducing the total wind mass to ~10°’ Mo. Although the 
signal peaks earlier it is still too late to explain our observations and the resultant spectrum is too 
high energy (peaking in the extreme UV; Extended Data Figure 6). We therefore conclude that an 
engine-driven kilonova is unlikely to explain the observed properties, whereas jet-ejecta 
interactions remain a plausible solution to enhance the early emission. 

In summary, we find that a compact binary merger would naturally account for most of the 
observed features of GRB211211A, from the onset of its kilonova to its environment and high- 
energy properties. The main challenge to this model remains the long duration of the prompt 


gamma-ray emission, requiring an active central engine for up to ~100 s. 
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Extended Data Table 1 — UV, optical and nIR observations of GRB211211A. Upper limits 


(u.l.) are 3 o. Values are corrected for Galactic extinction in the direction of the burst, Av = 0.047 
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mag’’. 
T - To Exposure Telescope Filter Magnitude Error 
(d) (s) AB (68% c.l.) 
0.0009 9 UVOT v 17:3 ul. 
0.0012 36 UVOT White 20.33 0.25 
0.0020 96 UVOT White 212 ul. 
0.042 195 UVOT u 19.75 0.13 
0.044 193 UVOT b 19.62 0.2 
0.046 185 UVOT White 19.6 0.06 
0.049 196 UVOT w2 19.59 0.12 
0.051 195 UVOT v 19:27 0.28 
0.053 196 UVOT m2 19.58 0.18 
0.056 196 UVOT wl 19.44 0.12 
0.058 78 UVOT u 19.38 0.16 
0.19 877 UVOT u 19575 0.07 
0.20 570 UVOT b 19.79 0.17 
0.24 3300 MITSUME g 19.85 0.15 
0.24 3300 MITSUME Re 20.3 0.19 
0.24 3300 MITSUME I, 20.2 ul. 
0.25 467 UVOT m2 20.5 0.16 
0.29 3300 MITSUME g 20.19 0.16 
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